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An  electro-thermal  cycle  life  model  is  developed  by  incorporating  the  dominant  capacity  fading  mech¬ 
anism  to  account  for  the  capacity  fading  effect  on  the  lithium  ion  battery  performance.  This  model  is 
comprehensively  validated  in  three  different  aspects:  (1)  Electrochemical  performance  with  different 
discharge  C_rates  (C/10, 1C,  3C)  and  different  ambient  temperature  (0  °C,  25  °C,  45  °C,  60  °C);  (2)  battery 
surface  temperature  and  center  temperature  during  6  Crate  constant  current  charge-constant  voltage 
charge-constant  current  discharge  (CC-CV-CC);  (3)  the  amount  of  capacity  fade  indicated  by  the 
discharge  potential  curves  of  different  cycle  conditions.  The  feasibility  of  using  the  dominant  factor  only 
to  represent  the  overall  capacity  fading  in  model  is  verified,  and  the  model  is  used  to  predict  OCV— SOC 
curves  of  cycled  battery.  To  recover  the  capacity  fade  of  cycled  battery,  either  charging  the  negative 
electrode  of  fully  discharged  battery  or  charging  the  positive  electrode  of  fully  charged  battery,  could  be 
an  effective  method,  and  possible  capacity  recovering  design  is  presented.  The  decreasing  trend  of 
thermodynamic  parameter  Ea  with  increasing  C_rate  in  the  capacity  fading  function  is  deduced  with 
increasing  battery  temperature  when  increasing  discharge  Crate. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Growing  concerns  have  been  paid  on  the  development  of 
lithium  ion  battery  and  its  application  on  electrical  vehicle  (EV) 
and  hybrid  electrical  vehicle  (HEV).  Manufacturer  of  both  lithium 
ion  battery  and  EV  need  a  better  understanding  of  lithium  ion 
battery  capacity  fading  mechanism,  and  have  a  major  interest  in 
predicting  and  improving  the  cycle  life  of  lithium  ion  battery  [1]. 
Many  Battery  management  systems  use  battery  open  circuit 
voltage  (OCV,  can  be  obtained  by  discharging  the  battery  at 
significantly  low  current)  to  predict  the  state  of  charge  (SOC)  of 
battery  [2,3],  however  OCV-SOC  curve  of  battery  changes  with 
battery  capacity  fading,  leading  to  inaccurate  estimation  of  SOC  for 
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cycled  battery  if  the  estimation  is  based  on  the  original  OCV-SOC 
curve.  So  the  prediction  of  OCV-SOC  curve  of  cycled  battery  is 
necessary  for  more  accurate  SOC  estimation.  The  capacity  fading  of 
lithium  ion  battery  is  caused  by  several  different  mechanisms 
associate  with  side  reactions,  leading  to  electrolyte  decomposi¬ 
tion,  passive  film  formation,  active  material  dissolution,  and  other 
phenomena  [4].  Recent  investigations  on  lithium  iron  phosphate 
battery  [5]  reveals  that  battery  capacity  is  affected  by  the  battery 
temperature,  depth  of  discharge  (DOD)  and  operating  current 
density. 

In  order  to  verify  capacity  fading  mechanisms  and  predict 
capacity  fading  within  battery,  capacity  fading  models  (Electro¬ 
chemical  model  [4],  Empirical  correlations  [6];  Equivalent  circuit 
model  [7])  have  been  presented  by  researchers.  We  are  interest  in 
electrochemical  models  and  empirical  correlations  in  this  paper. 
Electrochemical  models  incorporate  capacity  fading  governing 
equations  to  account  for  the  capacity  fading  effect  of  battery  using 
first  principle  model  introduced  by  Doyle  et  al.  [8].  Arora  et  al.  [4] 
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incorporated  multiple  capacity  fading  mechanisms  into  battery 
electrochemical  model,  and  suggested  another  model  [9]  consid¬ 
ering  the  lithium  deposition  overcharge  reaction  only  in  lithium 
ion  battery  with  carbon-based  negative  electrode,  the  growth  of 
SEI  layer  is  implemented  into  an  electrochemical  model.  SEI  layer 
has  been  widely  adopted  by  researchers  [10,11].  Empirical  corre¬ 
lations  are  developed  by  fitting  capacity  fading  rate  as  a  function 
of  factors  that  affect  the  capacity  fading  based  on  a  large  amount  of 
experiment  data.  Wang  et  al.  [6]  developed  a  cycle  life  model  by 
fitting  capacity  fading  as  a  function  of  battery  temperature, 
current  density  and  DOD  based  on  comprehensive  experimental 
data  [5]. 

Electrochemical  model  considering  all  the  capacity  fading 
mechanisms  has  been  proved  to  be  efficient  in  validating  capacity 
fading  mechanisms  as  discussed  above,  but  for  predicting  the  rate 
of  capacity  loss,  it  turns  out  to  be  tedious  work  and  it  is  difficult 
to  find  out  the  proportion  that  each  mechanism  contributes  to 
the  battery  capacity  loss  and  to  keep  a  balance  among  them  [12]. 
Therefore,  when  the  focus  is  on  predicting  the  rate  of  capacity 
loss,  implementing  the  dominant  capacity  fading  mechanism  into 
electrochemical  model  would  be  more  efficient  than  incorpo¬ 
rating  all  the  mechanisms.  Empirical  correlations  [6]  is  straight¬ 
forward  in  predicting  capacity  fading  rate  by  fitting  capacity 
fading  rate  as  a  function  of  battery  operating  conditions  (such  as 
temperature,  current  and  DOD),  however,  such  correlations  are 
unable  to  predict  corresponding  battery  discharging  performance 
and  OCV  curve  (which  is  essential  for  SOC  prediction).  In  this 
work,  we  try  to  predict  battery  cycle  life  and  battery  discharging 
performance  as  well  as  OCV  in  one  model  by  taking  advantages  of 
both  electrochemical  model  and  empirical  correlations  of  battery 
cycle  life.  In  order  to  take  into  account  of  the  temperature  effect 
on  capacity  fading,  an  electro-thermal  model  is  developed  than 
electrochemical  model.  Capacity  fading  mechanism  is  simplified 
that  the  dominant  mechanism  (active  lithium  loss)  is  chose  to 
account  for  the  overall  capacity  fading  process,  and  experimental 
results  [5]  and  correlations  [6]  are  incorporated  into  the  electro¬ 
thermal  model,  forming  an  electro-thermal  cycle  life  model  for 
lithium  ion  battery.  The  electro-thermal  performance  of  this 
model  is  validated  by  experimental  data  from  literature  [13,14], 
and  the  cycle  life  prediction  is  validated  by  experimental  data 
from  Wang  et  al.  [6].  The  feasibility  of  using  the  dominant 
mechanism  to  account  for  the  overall  capacity  fading  process  is 
verified  and  the  effect  of  operating  temperature  on  capacity 


fading  rate  is  evaluated  by  the  model.  The  model  can  be  used  for 
predicting  discharging  performance  and  OCV  of  cycled  battery, 
and  potential  design  consideration  and  capacity  renovation 
methods  are  suggested  and  discussed. 

2.  Model  development 

2.1.  Electrochemical  part 

The  electrochemical  model  for  LiFeP04  lithium  ion  battery  in 
this  paper  is  developed  based  on  the  pseudo  two-dimensional 
model  of  Doyle  et  al.  [8,15,16],  and  has  been  presented  in  our 
previous  work  [16].  Fig.  1  shows  a  schematic  calculation  domain  of 
one-dimensional  (ID)  battery  model  [8,15,17].  Two  inner  bound¬ 
aries  (Anode/separator  interface  boundaries  2  and  cathode/sepa¬ 
rator  interface  boundaries  3)  and  two  external  boundaries  (anode/ 
current  collector  interface  boundary  1  and  cathode/current 
collector  interface  boundary  4)  are  shown  in  the  figure.  The  spirally 
wound  cylindrical  battery  is  considered  as  a  lumped  single  cell,  so 
electrochemical  performance  is  assumed  to  be  homogeneous 
within  the  cylindrical  battery.  Average  battery  temperature  lavg  in 
the  2D  thermal  model  is  feedback  to  the  electrochemical  model  to 
account  for  the  temperature  effect  on  electrochemical  performance 
of  battery,  heat  generation  rate  Qavg  is  in  return  a  feedback  from 
electrochemical  model  to  thermal  model  for  calculating  the  battery 
temperature. 

2.1.1.  Electronic  charge  balance 

2. 1.1.1.  (a).  Solid  phase.  Electronic  charge  balance  for  solid  phase 
can  be  expressed  as 

V-(  — jcfv^)  =  — Sa>ijl0Cii  (1) 


where,  </>i  is  arbitrarily  set  at  zero  at  boundary  1 ;  At  boundary  4,  the 
charge  flux  is  set  to  be  equal  to  the  local  current  density  /;  At 
boundary  2  and  boundary  3  are  set  as  isolation,  since  there  is  no 
charge  flux,  Eq.  (4). 


Fig.  1.  Schematic  of  electro-thermal  model  geometry. 
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<Pi\x=o  =  0;  =  -J  (3) 

I X — Ln  +LS  +Lp 

-Kfvfa  I  =0;  -xfvft  1=0  (4) 

I X — Lf]  I  I X — Ln 


2A.1.2.  (b ).  Solution  phase.  The  electronic  charge  balance  governing 
equation  in  solution  phase  is  expressed  as: 
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Liquid-junction  potential  is  introduced  in  Eq.  (5)  with 
expression 
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where  v  is  the  thermodynamic  factor  relating  to  electrolyte  activity, 
it  is  temperature  and  concentration  dependent. 

<fc  is  taken  to  be  continuous  at  inner  boundaries  (boundary  2 
and  boundary  3),  while  it  is  set  to  be  isolate  boundary  at  external 
boundaries  (boundary  1  and  boundary  2)  since  there  is  no  flux  at 
these  boundaries,  Eq.  (8). 
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2.1.2.  Mass  balance 

2 A2A.  (a).  Solid  phase.  The  mass  balance  of  lithium  ions  in  an 
intercalation  particle  of  electrode  active  material  is  described  by 
Fields  law.y  is  the  distance  from  the  center  of  solid  sphere,  the  mass 
transport  within  solid  phase  can  be  described  as: 


dci  0 

-2--\ - n 

dt  0y2 


m  o 


(9) 


where  rp  is  the  particle  radius.  The  lithium  concentration  at  the 
surface  of  the  particles,  y  =  rp,  is  coupled  to  the  concentration  and 
flux  in  the  1 D  model  for  the  charge  and  material  transport  in  the 
electrolyte.  The  flux  at  the  center  of  sphere,  y  =  0,  is  set  to  zero, 
because  there  is  no  species  source. 


2.I.2.2.  (b).  Solution  phase.  Solution  phase  material  balance  for  LiF6 
dissolved  in  the  liquid  phase 

£2 :^  +  v{— DfVC2}  «4^(1-t+)  (10) 


Df  =  D2e ? 


(ID 


where,  jj0c,i  is  driven  by  overpotential,  rju  defined  as  the  difference 
between  solid  and  electrolyte  phase  potentials  minus  U,  the  ther¬ 
modynamic  equilibrium  potential  of  the  solid  phase. 

Vi  —  01  ,i  02, i  —  ^ref,i  0^) 

where  the  equilibrium  potential,  U,  is  taken  to  be  a  function  of  the 
solid  phase  lithium  ion  concentration  at  the  particle  surface.  The 
temperature-dependent  open  circuit  potential  of  electrode  i  is 
approximated  by  Taylor’s  first  order  expansion  around  a  reference 
temperature: 


U;  -  Uref  i  +  (t  -  Tref)  ^  (14) 

where  Ure ^  is  the  open  circuit  potential  under  the  reference 
temperature. 

In  Eq.  (12),  exchange  current  density,  z’o.i,  acts  as  a  bridge  con¬ 
necting  concentrations  in  both  solid  phase  as  liquid  phase: 

*0,i  ~  ^0,iC2C  (cl,max,i  —  cl,surf,i)  Cysurf,i  0^) 

where  k o,i  is  a  reaction  rate,  considered  temperature  dependent  in 
this  paper.  aa  and  ac  are  the  anodic  and  cathodic  transfer  coeffi¬ 
cients,  respectively. 

The  outputs  of  the  model  are  the  cell  potential,  current  density 
distribution,  species  and  concentrations  distributions.  The  cell 
potential  is  derived  by  following  expression: 

E  =  01  lx=Ln+Ls+Lp  -0llx=O  (I6) 


2.2.  Thermal  model 

There  are  3  domains  in  the  2D  thermal  model  geometry  in  Fig.  1 : 
1.  Battery  core  filled  with  electrolyte;  2.  Battery  body  filled  with 
battery  material  (electrodes,  current  collectors  and  separator);  3. 
Metal  can.  There  is  heat  generated  in  the  battery  body,  but  not  in 
battery  core  and  metal  can.  The  volume  averaged  heat  generation 
rate  QaVg  in  battery  body,  as  shown  in  Fig.  l,is  the  sum  of  the  local 
heat  generation  from  electrochemical  model,  including  reaction 
heat  generation  Qrea,  ohmic  heat  generation  Qohm  and  active  heat 
generation  QaCt,  heat  generation  in  current  collector  is  neglected 
because  the  electrical  conductivities  of  current  collectors  are 
significantly  larger  than  that  of  electrode  materials  that  heat 
generated  in  current  collectors  is  much  smaller  than  in  electrodes. 

Active  heat  generation  is: 

Qact  —  Sa  i Jloc  i  (0i  i  —  02,i  —  ^i)  (17) 

Reaction  heat  generation  is: 


The  flux  of  liquid  species  is  set  to  zero  at  external  boundaries 
(boundary  1  and  4),  liquid  species  flux  and  species  concentration  Qrea  =  Sa ,i Jioc jT 
are  taken  to  be  continuous  at  inner  boundaries  (boundary  2  and 

boundary  3).  Ohmic  heat  generation  is: 


2 A. 3.  Electrochemical  kinetics 

The  local  current  per  active  material  area  is  calculated  using  the 
Bulter-Volmer  equation: 


Qohm  =  <^<£1  •  V0!  +  KeffV</>2  '  V<fo  + 

Ce 

So  the  volume  averaged  heat  generation  rate  is 


(18) 


(19) 
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Qavg  —  Qact  +  Qrea  +  Qohm  (20) 

The  energy  balance  of  the  2D  geometry  is  shown  as  follows: 
pC/§:+V-(-kVT)  =  ^ 

p,  Cp  and  k  are  the  volume  averaged  density,  heat  capacity  and 
thermal  conductivity,  respectively. 

According  to  Newton’s  cooling  law  and  radiation  law,  the 
boundary  condition  for  energy  balance  is  expressed  as: 

=  -h(Tamb  - T)  -  £(T lb  -  jA)  (21) 

where  h  is  lumped  heat  transfer  coefficient,  Tam b  is  the  ambient 
temperature,  and  e  is  the  blackness  of  the  battery  surface  (e  =  0.8 
[18]). 

2.3.  Capacity  fading  part 

The  capacity  of  lithium  ion  battery  is  limited  by  the  amount  of 
active  lithium  [5].  The  instability  of  carbon  negative/electrolyte 
interface  is  suspected  to  be  the  source  of  active  lithium  loss. 
However,  there  is  slight  change  in  the  positive  electrode  [5].  This 
paper  focuses  on  predicting  discharging  performance  and  OCV  of 
cycled  battery  in  a  simple  way,  so  we  incorporate  only  the  domi¬ 
nant  mechanism  (active  lithium  loss)  into  the  model  for  cycle  life 
prediction,  the  feasibility  of  such  simplification  will  be  discussed 
later.  Therefore,  in  this  model,  the  change  of  amount  of  active 
lithium  (dominant  mechanism)  is  expressed  by  varying  the  initial 
state  of  charge  (SOCo,n)  of  negative  electrode  to  account  for  the 
battery  capacity  fading  process. 

As  mentioned  above,  Wang  et  al.  [6]  provided  a  function  of 
capacity  fading  rate  of  26,650  cylindrical  lithium  iron  phosphate 
during  C/2  discharging  at  different  ambient  temperature. 

Qioss  =  30330exp(:g^A°552  (  22) 

In  this  work,  the  above  function  is  integrated  into  SOCo,n  to 
account  for  the  capacity  fading  effect  in  the  electro-thermal  model. 

SOConN  =  S0Co,n  (l  -  30330exp  (' g^”)^'552)  (23) 

where  Ah  is  expressed  as 

Ah  =  (cycle_number)  x  (DOD)  x  (fulLcelLcapacity)  (24) 

By  imputing  the  cycle  number  N,  this  model  would  be  able  to 
calculate  the  capacity  loss  of  battery. 

3.  Model  parameters  and  model  validation 

3 A.  Model  calibration 
3.1.1.  Model  parameters 

The  calculations  of  model  were  performed  using  the  finite 
element  commercial  software  COMSOL  MULTIPHYSICS®  (Version 
3.5a).  Model  parameters  are  from  manufacturer’s  data,  literature 
and  estimation.  Parameters  for  a  2.3  Ah  LiFePCH  battery  are  listed  in 
Tables.  1  and  2,  including  design  specifications,  lithium  ion 
concentration  parameters,  Kinetic  and  transport  properties  and 
thermal  properties.  Kinetic  reaction  rate  constant  at  reference 
temperature  of  cathode  (/<o,p)  in  Table  1  is  estimated  by  matching 
the  simulated  discharge  potential  with  experiment  data.  Specific 
heat  capacity  and  density  of  the  battery  body  in  Table  2  are  the 


Table  1 

Model  parameter  for  a  2.3  Ah  LiFeP04  cylindrical  26650  type  battery  electrochemical 
sub-models. 


Anode 

Separator 

Cathode 

Design  specifications  (geometry  and  volume  fractions) 

Aceii  (m2) 

0.16943 

£l,i 

0.55a 

- 

0.43a 

£2,i 

0.33a 

0.54a 

0.332a 

U  (pm) 

34a 

30a 

70a 

Ri(pm) 

0.0365a 

- 

3.5a 

Lithium  ion  concentrations 

C2o  (mol  m-3) 

1200a 

Cmax.i  (mol  m-3) 

31370b 

- 

22806b 

SOC0,i 

0.8b 

- 

0.03b 

Kinetic  and  transport  properties 

aa.i,  «c.i 

0.5 

- 

0.5 

7i 

1.5 

1.5 

1.5 

D2  (m2  s-1) 

Eq.  (25) 

Di0,i  (m2  s-1) 

3.9  x  10-14c 

- 

1.18  x  10-18b 

£a,D,i  (kj  mol"1) 

35d 

- 

35[b] 

/<o,i  (m2  5  mol-0-5  s-1) 

3  x  10-lle 

- 

1.4xl0-12f 

£a,k,i  (kj  mor1) 

20g 

- 

30b 

K\  (S  nr1) 

100c 

- 

0.5h 

k2  (S  m-1) 

Eq.  (26) 

t+ 

0.363c 

V 

Eq.  (27) 

Tref(I<) 

298.15 

F  (C  mol-1) 

96,487 

a  Manufacturer. 
b  Ref.  [13]. 
c  Ref.  [19]. 
d  Ref.  [20]. 
e  Ref.  [21]. 
f  Estimate. 
g  Ref.  [22]. 
h  Ref.  [23]. 


volume  averaged  value  of  electrode  materials,  current  collectors 
and  separator.  Bulk  thermal  conductivity  is  experimentally 
measured  to  be  0.2  W(m  K)-1  rather  than  using  the  volume  aver¬ 
aged  value  because  there  are  contact  resistances  which  are  difficult 
to  measure  between  two  different  layer  of  battery  materials. 

3.1.2.  Temperature  dependent  parameters 

The  data  of  lithium  ion  diffusivity  in  solution  phase  D2,  electric 
conductivity  in  solution  phase  k2  and  thermodynamic  factor  v  are 
adopted  from  published  literature  Ref.  [28].  The  solvent  mixture  for 
the  commercial  battery  tested  is  unknown,  but  we  assume  that  Eqs. 
(25)— (27)  are  applicable  in  this  commercial  battery  and  use  these 
data  in  our  calculation  even  though  the  solvent  mixture  of  the 
battery  might  be  different  from  that  in  [28]: 


Table  2 

Thermal  properties  for  a  2.3  Ah  LiFeP04  cylindrical  26650  type  battery. 


Thermal  properties 

Id  (W(m  I<)-1) 

Pi  (kg  m  3) 

CpilpgK)-1) 

Anode 

1.04a 

2500b 

800b 

Separator 

la 

1200b 

800b 

Cathode 

1.48a 

1500b 

800b 

Cu 

401 

8900 

3440 

Al 

237 

2700 

2420 

Battery  material 

0.2C 

2101 

1014 

Volume  averaged 
Electrolyte 

0.6d 

1130d 

2055d 

Can 

14e 

7500e 

460e 

a  Ref.  [24]. 
b  Ref.  [25]. 
c  Measure. 
d  Ref.  [26]. 
e  Ref.  [27]. 
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=  1  XlO  410  443  r-22950-0.005c  22  Xl°  4°  (25) 

k2  =  1  x  10_4c(  -  10.5  +  0.074T  -  6.69  x  10“5T2 

+6.68  x  10-4c-  1.78  x  10~5c7 
+2.8  x  10_8cT2  +  4.94  x  lO^c2 

-8.86  x  W~w(?T)2  (26) 

v  =0.601  -0.24Vl0-3c 

r  , - 1  (27) 

+  0.982  1  -  0.0052(7  -  294.0)  V  10-9c3 

Diffusivity  coefficient  and  reaction  rate  constants  follow  the 
following  Arrhenius  formula  [29]: 


SOC 

Fig.  2.  Open  circuit  potential  of  positive  electrode  (Up)  and  negative  electrode  (Un). 
Lithium  iron  phosphate  entropy  is  also  from  Ref.  [23]  : 


A  = 


A),iexP 


D1>n  =  3.9  x  10 
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V  R  1 

(298.15 

T 

(28)  =  ~  °-35376*8  +  1  -3902x7  +  1.9635X5  -  0.9871  6X4 

+  0.28857X3  -  0.046272X2  +  0.0032158x 
-  1.9186  x  105 

(29)  (36) 


Di  p  =  1.18  x  10_18exp^ 
/<l  n  =  3  x  10_11exp( : 


/35000 
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l  R  1 
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( 20000 

(  1 

V  R  1 

(298.15 
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3.2.  Model  validation 

(30) 

This  electro-thermal  cycle  life  model  is  validated  from  electro¬ 
chemical  performance,  thermal  performance  and  cycle  life 
perspective.  Experimental  data  are  from  different  experiment  done 

(31)  by  different  researchers  [6,13,14]  with  the  same  type  of  battery 
(26650C  lithium  iron  phosphate  battery,  2.3  Ah). 


fei,P 


1.4  x  10  12exp(-3y)exp 


/30000 

(  1 

!Y> 

l  R  1 

(298.15 

T 

(32) 


In  order  to  match  the  simulation  result  with  experimental  data, 
we  add  a  ratio  of  exp(-3y)  into  the  formula  of  reaction  rate  of 
positive  electrode.  This  ratio  is  specific  to  our  model,  and  does  not 
necessarily  apply  to  other  models. 


3.1.3.  Electrode  thermodynamic  properties 

State  of  charge  (SOC)  of  LixC6  negative  electrode  is: 

SOCn  =  x  =  n/C|  max  n  (33) 

And  state  of  charge  of  LiyFeP04  positive  electrode  is: 


3.2.1.  Electrochemical  performance 

The  electro-thermal  cycle  life  model  is  first  validated  from  the 
electrochemical  performance  perspective.  Experiment  data  of 
battery  potential  against  discharge  capacity  at  different  ambient 
temperature  (0  °C  [6],  25  °C  [13],  45  °C  [6,13],  60  °C  [6])  with  C/2 
current  are  compared  with  the  simulation  results.  The  good 
agreements  in  Fig.  3  shows  that  this  electro-thermal  cycle  life 
model  is  able  to  simulate  battery  discharge  performance  at 
different  ambient  temperature.  To  further  validate  the  electro¬ 
chemical  performance  during  charge  and  discharge  with  different 
current  density,  the  simulation  results  are  compared  with  the 
experiment  data  [13]  of  battery  charge  and  discharge  with  different 
current  density  (C/10,  1/2C,  1  C)  at  45  °C  in  Fig.  4  and  Fig.  5, 
respectively. 


SOCp  —  y  —  C|p/C|maXp  (34) 

The  open  circuit  potential  of  negative  electrode  Un  and  positive 
electrode  LTp  are  from  Ref.  [13],  Fig.  2.  Point  B  is  the  characteristic 
point  of  open  circuit  potential  curve  of  negative  electrode  (l/n),  DE  is 
the  effective  SOC  range  of  positive  electrode,  and  AC  is  the  effective 
SOC  range  of  negative  electrode. 

Lithium  graphite  entropy  is  from  Ref.  [23]: 


3.2.2.  Thermal  characterization 

Forgez  et  al.  [14]  experimentally  measured  the  battery  surface 
temperature  and  center  temperature  during  6  C  constant  current 
charge-constant  voltage  charge-constant  current  discharge 
process  (CC-CV-CC).  To  evaluate  the  thermal  performance  of  the 
electro-thermal  cycle  life  model,  the  current  profile  of  Forgez  et  al. 
[14]  is  implemented  into  the  model  for  simulation.  The  comparison 
of  simulation  results  with  experiment  data  is  shown  in  Fig.  6.  The 
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344.1347148 


exp(-32.9633287x  +  8.31671 1484) 

1  +  749.0756003exp(-34.7909964x  +  8.887143624) 


-  0.8520278805X  +  0.36229929x2  +  0.2698001 697 
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Fig.  3.  0  °C  [6],  25  °C  [13],  45  °C  [6,13],  60  °C  [6],  C/2  discharge  validation. 


experimental  data  of  temperature  include  the  center  temperature 
of  battery  measured  by  a  thermal  couple  inserted  into  the  center  of 
the  battery,  and  the  surface  temperature  measured  by  a  thermal 
couple  stuck  to  the  battery  surface  [14].  Temperature  decreases 
during  the  constant  voltage  charging  period,  probably  because 
there  is  lower  heat  generation  rate  during  that  period  than  constant 
current  charging  period  and  constant  current  discharging  period. 
The  simulation  results  agree  well  with  experimental  data  during 
the  CC-CV— CC  process,  showing  that  this  model  is  capable  of 
predicting  battery  temperature. 

3.2.3.  Capacity  fading  validation 

Finally,  the  simulation  results  of  discharge  performance  at 
different  ambient  temperature  and  different  cycle  number  are 
validated  against  the  experiment  data  from  Wang  [6],  Fig.  7  and 
Fig.  8.  It  has  to  be  mentioned  that  the  cycle  life  function  imple¬ 
mented  in  this  electro-thermal  model  is  also  adopted  from  the 
work  of  Wang  [6].  Therefore  the  accuracy  of  capacity  prediction 
in  this  model  (related  to  the  end  point  correspond  to  the 
discharge  end)  would  be  the  same  as  the  function  of  Wang’s. 
However,  the  good  agreement  of  discharge  curve  shapes  is 
ensured  by  the  implementing  method  of  capacity  fade  function 
into  the  model,  this  method  would  be  further  discussed  in  the 
following  section. 


4.  Results  and  discussion 

4.1.  Assumption  validation  and  OCV  prediction 

As  we  assume  that  the  dominant  factor  of  capacity  fading  can 
represent  the  overall  capacity  fading  phenomenon,  and  we  use  the 
variation  of  SOCo,n  alone  to  simulate  the  capacity  fading  in  the 
model  for  the  sake  of  simplicity,  though  the  discharging  perfor¬ 
mances  of  cycled  battery  have  been  validated  in  the  previous 
section,  we  still  need  to  verify  the  feasibility  of  making  the 
assumption  more  clearly.  As  shown  in  Fig.  9,  the  simulation  result 
of  discharging  potential  of  0.5  C  cycled  272  times  at  45  °C  matches 
well  with  the  experiment  data.  The  characteristic  point  A  nearly 
overlap  between  simulation  result  and  experiment  data,  while 
there  is  a  slight  discrepancy  of  the  position  of  characteristic  point  B 
between  simulation  result  and  experiment  data.  The  discrepancy  of 
B  may  be  due  to  our  assumption  of  using  the  dominant  factor  to 
represent  the  overall  capacity  fading  phenomenon.  However,  as  far 
as  discharging  potential  is  concerned,  the  small  discrepancy  in 
discharging  potential  is  acceptable. 


Fig.  6.  Surface  temperature  and  center  temperature  of  battery  during  6C  (CC-CV-CC) 
charge  and  discharge  [10]. 
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Fig.  7.  Validation  of  capacity  fading  at  45  °C  [6]. 


Fig.  9.  45  °C,  1/2  Cirate  discharge,  cycle  number  272,  verification  of  discharging 
potential. 


As  introduced  before,  OCV-SOC  curve  in  BMS  should  be  upda¬ 
ted  than  the  original  one  to  ensure  accuracy  in  SOC  estimation.  As 
OCV-SOC  curve  can  be  experimentally  obtained  by  discharging  the 
battery  at  low  current,  we  try  to  simulate  the  OCV-SOC  curves  of 
cycled  battery  by  discharging  the  battery  with  1  /50  C  current  in  the 
model.  Simulation  results  are  shown  in  Fig.  10,  the  OCV  curves  of 
different  cycle  numbers  differ  a  lot,  which  in  return  verify  the 
significance  of  using  updated  OCV  curves  than  the  original  one  in 
BMS  for  SOC  estimation. 

4.2.  Capacity  recovering 

As  shown  in  Fig.  11,  due  to  the  loss  of  active  lithium  in  cycled 
battery,  for  a  fully  discharged  battery,  the  SOC  of  the  positive 
electrode  moves  from  E  (the  original  position)  to  Ei,  and  for  a  fully 
charged  battery,  the  SOC  of  the  negative  electrode  moves  from  C 
(the  original  position)  to  Cp 

To  recover  the  capacity  in  cycled  battery,  two  methods  can  be 
utilized.  When  the  cycled  battery  is  fully  charged  (point  Ci  in 
negative  electrode  and  point  D  in  positive  electrode  in  Fig.  11), 
further  charging  is  impossible  since  there  is  not  further  supply  of 
lithium  ions  from  positive  electrode  to  negative  electrode.  There¬ 
fore,  the  first  method  is  to  charge  the  positive  electrode  alone  from 
D  to  D2,  this  charging  process  can  be  conduct  by  charging  the 
positive  electrode  with  a  lithium  metal  electrode  (lithium  metal 
electrode  acts  as  lithium  ion  source  for  positive  electrode).  After  the 
process,  the  battery  can  be  further  charged,  since  there  are  more 


lithium  ions  in  positive  electrode  to  transport  to  the  negative 
electrode.  The  final  condition  of  fully  charged  battery  is  point  D  in 
positive  electrode  and  point  C2  in  negative  electrode.  And  fully 
discharged  battery  is  point  E2  in  positive  electrode  and  point  A  in 
negative  electrode.  Capacity  before  recovering  is  ACi,  and  after 
recovering  is  AC2,  capacity  of  C1C2  is  recovered  to  the  battery.  The 
schematic  setup  of  the  recovering  process  is  shown  in  Fig.  12, 
a  lithium  metal  electrode  is  connected  with  the  positive  electrode 
current  collector  of  the  disclosed  battery  by  the  external  power 
supply  and  the  same  electrolyte  as  the  battery.  Lithium  metal 
serves  as  the  lithium  ion  sources  and  the  electrolyte  acts  as 
a  channel  for  lithium  ions  to  transport  to  and  insert  into  the  positive 
electrodes.  In  this  process,  additional  active  lithium  is  supplied  to 
the  cycle  battery. 

On  the  other  hand,  when  the  battery  is  fully  discharged  (nega¬ 
tive  electrode  at  point  A  and  positive  electrode  at  Ei ),  the  capacity  is 
limited  by  the  negative  electrode  that  its  SOC  is  0  (point  A)  and 
cannot  further  provide  lithium  ions  to  insert  into  the  positive 
electrode,  Fig.  11.  Therefore,  another  method  for  recovering  battery 
capacity  is  to  charge  the  negative  electrode  alone  to  move  point  A 
rightward  to  A3,  this  process  can  be  conduct  by  charging  the 
negative  alone  with  a  lithium  metal  reference  electrode.  After  this 
process,  the  battery  can  be  further  discharged,  since  there  are  more 
lithium  ions  in  negative  electrode  to  provide  to  insert  into  the 
positive  electrode.  The  final  condition  of  fully  discharged  battery 
can  be  point  E3  in  positive  electrode  and  point  A  in  negative 
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Fig.  11.  Schematic  of  battery  capacity  fading  and  recovering. 

electrode,  and  fully  charged  battery  can  be  point  D  in  positive 
electrode  and  point  C3  in  negative  electrode.  Capacity  before 
recovering  is  ACi,  and  after  recovering  is  AC3,  capacity  of  C1C3  is 
recovered  to  the  battery.  The  schematic  setup  of  the  recovering 
process  is  shown  in  Fig.  13. 

Considering  the  safety  issue  in  the  recovering  process  of 
capacity  fade,  charging  the  negative  electrode  of  fully  discharged 
battery  would  be  a  better  choice  than  charging  the  positive  elec¬ 
trode  of  fully  charged  battery,  because  short  circuit  or  other  hazards 
would  cause  more  serious  safety  problem  for  fully  charged  battery 
than  fully  discharged  battery.  As  shown  in  Fig.  14(a),  the  core  of 
a  cylindrical  battery  is  filled  with  electrolyte,  a  lithium  metal  rode 
electrode  can  be  inserted  into  the  core,  it  would  be  possible  espe¬ 
cially  for  battery  with  larger  diameter.  Charging  the  negative 
electrode  of  a  fully  discharged  battery  may  be  realized  by  the  setup 
sketched  in  Fig.  14(b). 

4.3.  Self-heating  and  temperature  effect 

Capacity  fading  rate  is  a  function  of  temperature,  while  the 
battery  temperature  would  be  rose  by  the  self-heating  effect  of 


I  u-  1 

Electrolyte 


Fig.  12.  Schematic  setup  of  recovering  battery  capacity  by  charging  positive  electrodes. 


7N — 

Li" 


Electrolyte 


Fig.  13.  Schematic  setup  of  recovering  battery  capacity  by  charging  negative 
electrodes. 

battery,  leading  to  higher  temperature  of  battery  and  higher 
capacity  fading  rate.  Capacity  fading  rate  can  be  expressed  as  [8]: 

Qioss  =  Bexp(^)(ah)z  (38) 

where  the  pre-exponent  factor  B  decreases  with  increasing  C_rate, 
but  is  difficult  to  quantitatively  describe  this  relationship  by  using 
a  simple  mathematical  correlation;  The  index  number  z  is  almost 
constant  to  0.55;  while  the  Ea  decrease  with  increasing  C_rate,  and 
is  fitted  into  a  function  of  C_rate  as  follows  [6]: 

Ea  =  -31700  +  370.3  x  C_rate  (39) 

Flowever,  as  a  thermodynamic  term  for  battery,  Ea  is  supposed 
to  be  a  constant  rather  than  a  function  of  C_rate.  The  decreasing 
trend  of  Ea  with  C_rate  in  Eq.  (38)  is  deduced  to  the  lower  T  value 
used  in  the  function.  Wang  et  al.  [6]  used  the  battery  surface 
temperature  in  the  Eq.  (38),  while  the  surface  temperature  may  not 
be  representative  for  the  battery.  Due  to  the  high  heat  transfer 
coefficient  around  the  battery  surface  in  experiment,  the  surface 
temperature  of  the  battery  is  close  to  the  ambient  air  temperature, 
especially  for  those  conditions  in  which  the  ambient  air  tempera¬ 
ture  is  higher  [6],  such  as  25  °C  or  45  °C,  while  the  temperature 
within  battery  differ  from  the  surface  temperature  and  is  difficult  to 
measure. 

Simulation  results  (Fig.  15)  of  battery  discharged  with  different 
Clrate  at  ambient  temperature  of  45  °C,  and  heat  transfer  coeffi¬ 
cient  h  at  70  W  m~2  K-1  show  that  the  surface  temperature  of 
battery  is  generally  lower  than  the  average  temperature,  and  the 
difference  between  battery  surface  temperature  and  average 
temperature  increases  with  Clrate  and  discharge  capacity  as  shown 
in  Fig.  10.  Therefore,  battery  surface  temperature  (Tsurface)  under¬ 
estimates  the  battery  temperature.  To  better  capture  the  influence 
of  temperature  on  the  capacity  fading,  using  the  volume  average 
temperature  as  a  characteristic  temperature  would  be  a  better 
choice,  and  the  Ea  in  the  capacity  function  is  likely  to  be  a  constant 
at  different  discharging  CLrates.  In  sum,  surface  temperatures  used 
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current  collector  ^ 
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and  current  collector 


Electrolyte 


3  Top  view  of  disclosed  battery 


b  Side  view  of  disclosed  battery 


Fig.  14.  Schematic  setup  of  recovering  cylindrical  battery  capacity  by  charging  the  negative  electrodes. 


Fig.  15.  Temperature  difference  between  average  temperature  and  surface  tempera¬ 
ture,  Tavg  Tsurface- 

in  in  Wang’s  [6]  capacity  fading  function  underestimate  the 
temperature  effect  of  capacity  fading  rate.  However,  due  to  the 
measurement  limitation,  the  measurable  temperature  is  the 
surface  temperature,  it  is  a  compromise  to  use  the  surface 
temperature  as  the  characteristic  temperature  of  the  battery. 

5.  Conclusion 

An  electro-thermal  cycle  life  model  is  develop  by  implementing 
capacity  fading  effect  in  electro-thermal  model  of  cylindrical 
lithium  ion  battery,  this  model  is  able  to  simulate  the  discharging 
performance  during  different  discharge  cycles,  predicting  battery 
temperature,  as  well  as  predicting  capacity  loss  at  different  cycle 
number.  The  electro-thermal  cycle  model  is  comprehensively 
validated  in  the  aspects  and  the  simulation  results  well  agree  with 
experimental  data. 

It  is  feasible  to  use  the  dominant  factor  to  represent  the  overall 
capacity  fading  of  battery  for  simplicity,  and  the  accuracy  of 
simulation  results  of  discharging  potential  is  acceptable.  The  model 


can  be  used  to  predict  updated  OCV  curves  which  are  necessary  for 
accurate  estimation  of  SOC  for  cycled  battery  in  BMS. 

To  recover  the  capacity  loss  of  cycled  battery,  additional  lithium 
ion  should  be  provided  to  the  battery.  Either  charging  the  negative 
electrode  of  fully  discharged  battery  or  charging  the  positive  elec¬ 
trode  of  fully  charged  battery,  would  be  possible  method  to  achieve 
active  lithium  compensation.  Charging  the  negative  electrode  of 
fully  discharged  battery  would  be  a  better  choice  considering  the 
safety  issue.  Possible  capacity  recovering  setup  has  been  proposed. 

The  temperature  predicting  capability  is  validated  against  the 
experimental  data  of  battery  surface  temperature  and  center 
temperature  during  6  C_rate  constant  current  charge-constant 
voltage  charge-constant  current  discharge  (CC-CV-CC).  Temper¬ 
atures  predicted  by  the  model  for  different  C_rate  show  that  battery 
surface  temperature  underestimate  the  temperature  of  battery,  the 
decreasing  trend  of  Ea  with  increasing  C_rate  is  deduced  to  the 
increasing  battery  temperature  with  increasing  discharge  C_rate, 
and  half  radius  temperature  or  volume  average  temperature  would 
be  a  better  parameter  to  characterize  the  battery  temperature  and 
would  be  able  to  better  help  evaluate  temperature  effect  on 
capacity  fading  rate. 
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Nomenclature 

Aceii:  area  of  the  positive  (both  sides)  that  has  opposing  negative  (m2) 

C\X-  reduced-lithium  in  active  material  (mol  m-3) 

Cimax,i- ’  maximum  concentration  (mol  m-3) 

CPii:  heat  capacity  (J(kg  K)-1) 

Djx  solid  phase  diffusivity  (m2  s-1) 

DjoX  s°lid  phase  diffusivity  at  reference  temperature  (m2  s-1) 

Eq.dX  active  energy  for  diffusion  (kj  mol-1) 

Ea,k,i-  reaction  active  energy  (kj  mol-1) 
h:  heat  transfer  coefficient  (W  m-2  K-1) 

I:  cell  current  density  based  on  Aceii  (A  m-2) 
joX  exchange  current  density  (A  m-2) 
jl0Cii:  local  current  density  (A  m-2) 
koj:  reaction  rate  constant  (m2  5  mol-0  5  s-1) 
kj.'  thermal  conductivity  (W(m  I<)-1) 

U:  thickness  (pm) 

N:  cycle  number 

Qact'.  active  heat  generation  (W  m-3) 

Qohm :  ohmic  heat  generation  (W  m-3) 


Qrea:  reaction  heat  generation  (W  m-3) 

Qioss •  ratio  of  capacity  loss  over  full  cell  capacity 

r:  radius  distance  variable  of  particle  (m) 

rp:  characteristic  particle  radius  of  electrode  particles  (pm) 

SaX  specific  surface  area  (m-1) 

SOCoX  initial  state  of  charge 
t:  time  (s) 

t+:  transference  number  of  Li  ion  species  dissolved  in  liquid 
T:  absolute  temperature  (K) 

TamX  ambient  temperature  (K) 

U:  thermodynamic,  open  circuit  voltage  (V) 
v:  thermodynamic  factor  relating  to  electrolyte  activity 
x:  distance  variable  through  a  cell  component  (m) 
y:  dimensionless  radial  distance  of  particles 
x:  average  composition  variable  of  Li  in  Li*C6 
y:  average  composition  variable  of  Li  in  LiyFeP04 

Greek  letters 

aax  transfer  coefficient  for  anodic  current 

acx  transfer  coefficient  for  cathodic  current 

e:  blackness  of  battery  surface 

£jX  active  material  volume  fraction 

£2X  volume  fraction 

0i:  electric  potential  (V) 

jii  Bruggeman  tortuosity  exponent 

k:  ionic  or  electronic  conductivity  (S  m-1) 

Pi :  density  (kg  m-3) 

<Tj.‘  solid  phase  conductivity  (S  m-1) 

Subscripts,  superscripts  and  acronyms 
0:  initial  or  equilibrated  state 
amb:  ambient  (temperature) 
cell:  battery 
n:  negative  electrode 
p:  positive  electrode 
s:  separator 

ref:  reference  composition  of  relative  to  a  Li/lithium  ion  reference  electrode  or 
reference  temperature  for  Arrhenius  formula 


